Proposal for measurement of harmonic oscillator Berry phase in ion traps 
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We propose a scheme for measuring the Berry phase in the vibrational degree of freedom of a 
trapped ion. Starting from the ion in a vibrational coherent state we show how to reverse the sign of 
the coherent state amplitude by using a purely geometric phase. This can then be detected through 
the internal degrees of freedom of the ion. Our method can be applied to preparation of Schrodinger 
cat states. 
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When the Hamiltonian of a quantum system is varied 
adiabatically in a cyclic fashion, the state of the system 
acquires a geometrical phase in addition to the usual dy- 
namical phase. This effect, discovered by Berry (and 
generalized in various ways has been widely tested 
for 2-state systems || and attracted interest from a vari- 
ety of fields (HJI]. However, the Berry phase has not been 
experimentally measured for quantum harmonic oscilla- 
tors, though some theoretical calculations exist [|6|-^0|. 
The reason for this might be the fact that for the sim- 
plest case, namely for adiabatic displacement of an oscil- 
lator state in phase space, the Berry phase is indepen- 
dent of the state ||, and thereby undetectable. How- 
ever, when a squeezing Hamiltonian is switched on, and 
the squeezing parameter is varied, there would be a de- 
tectable Berry phase For an initial Fock state \n) 
which undergoes squeezing, the Berry phase after a cycle 
is — (n+ 1/2) times the classical Hannay angle H^]. The 
squeezed states of the electromagnetic field would have 
been a natural candidate to test this kind of phase, but 
they are not stable enough for an adiabatic evolution. 
However, the vibrational mode of a trapped ion has been 
a fertile ground for the preparation of long lived non- 
classical states of a harmonic oscillator In this 
letter, we derive a Berry phase formula for a certain adi- 
abatic evolution of a joint state of the internal levels of a 
trapped ion and its vibrational motion. Despite being the 
phase gained by a joint state, its value is fundamentally 
dependent on the harmonic oscillator nature of the vibra- 
tional mode. We propose a scheme to detect this phase 
which is feasible with current technology. It is worth- 
while to mention that based on calculations of Ref. ||] 
(which where tested in systems other than quantum har- 
monic oscillators ) , there has been an earlier attempt 
to detect the nonadiabatic geometric phase in ion traps 
by applying a set of four squeezes to the vibrational state 
p7[ . Here we propose a way of detecting the harmonic 
oscillator version of Berry's original adiabatic geometric 
phase. 

Consider the Hamiltonian 



H = H a + Hf, 



(1) 



g a e^\e){g\a + h.c. 



(2) 



Hb = 9b\e)(g\a^ + h.c. . 
In the above, |e) and \g) are two states of a qubit, 



(3) 
it 

and a are the creation and annihilation operators of a 
harmonic oscillator, g a and gb are unequal positive in- 
teraction strengths (say g a > gb) and <f> is an arbitrary 
phase factor. The motivation for choosing this Hamilto- 
nian is the possibility of its physical implementation and 
this will be described later. If the phase <f> is slowly varied 
over a complete loop (so that the adiabatic approxima- 
tion holds true), there will be a nontrivial Berry phase 
acquired by an eigenstate of the Hamiltonian H . We now 
proceed to calculate this. We transform the Hamiltonian 
as 



H = S{e)^HS{e) 



(4) 



where S{e)' i aS{e) = acosh(r) — sinh(r)e ie is a squeez- 
ing transformation with squeezing parameter e = re l6 . 
If we chose the squeezing strength r = tanh -1 gb/g a and 
the squeezing phase 8 = —<f>, the transformed Hamilto- 
nian will be 



H =il(\e){g\a + a*\g)(e\) 



(5) 



where Vt — g a cosh(r) — gb sinh(r). H is the well known 
resonant Jaynes-Cummings Hamiltonian [pi . A similar 
transformation to arrive at a Jaynes Cummings Hamil- 
tonian has been used in Refs. [|l9|,|o). The eigenstates of 
H are 



|*±> = -^(|3)|n + l>±|e)|n». 



(6) 



This implies that the eigenstates of our original Hamil- 
tonian H are 



%) 
v/2 



(| fl >|n + l)±|e)|n)). 



(7) 



where 



The states S(e)\n) featuring in the above expression are 
the squeezed number states [M . Now we can proceed to 



1 



calculate the Berry phase from the instantaneous eigen- 
states |**(e)) of H. 

The expression for the Berry phase for an adiabatic 
cyclic evolution of a Hamiltonian ff(Ft) in parameter 
space R is given in terms of the instantaneous eigenstates 
|n,R) of the Hamiltonian as 

7„ = ijdR(n,R\V R \n,R). (8) 

Using Eq.(^|) in the above equation we get 

ln =ijde{^\S(e)^ e S(e)\^) 

= ~Jde(n+l\S(e)tVcS(e)\n+l) 

+ l -jde{n\S{e)^ e S{e)\n). (9) 

If the modulus r of the parameter e is kept constant 
throughout the evolution, then using the expression for 
(n\S(e)^ e S(e)\n) from Ref. § we get 



7„ = — 2n(n + 1) sinh 2 r. 



(10) 



Note that the Berry phase 7„ is same for both the states 
|\I/+(e)) and ^"(e)). However, the dynamical phase 
/3jf is exactly opposite for the eigenstates |\l/+(e)) and 
|\I/~(e)) and is given by 



IK 



dt(*t\S(eyHS(e)\*t) 
T^Vn + IT, 



(11) 



where T is the time period of the cyclic variation of e. 
Thus one can make the dynamical phase completely van- 
ish after two cycles by changing the state |^+(e)) to 
or vice versa after one cycle. Under such cir- 
cumstances, the only contribution to the phase of the 
system will be geometrical. 

Let us now describe how the Hamiltonian H of Eq. ([!]) 
can be physically realized. Recently, ion tra ps h ave been 
a very active field of both experimental [jfl - 14 and the- 
oretical Jl5| , ^9pC| , |22] | research. Consider a single A three- 
level ion with two hyperfine ground states and an excited 
state (such as 9 Be + @) in a harmonic trap of frequency 
v . The two ground states, labelled as |e) and \g), are sep- 
arated in frequency by an amount uj which is much less 
than their separation from the excited state. The motion 
of the ion is modified by its interaction with two pairs 
of travelling- wave laser beams whose frequencies are de- 
tuned from the excited state. The first pair of lasers have 
frequencies uli and u>L2 which satisfy luli — ^L2 = lo + v 
and the second pair of lasers have frequencies uj^ and 
lul4 which satisfy u>l3 — = — v. We require that 
the pair of lasers LI and L2 differ in frequency from the 
pair L3 and LA by by an amount much larger than loq. 
We also assume that LI differs in phase from the rest of 



the beams by an amount (j> (this phase will need to be 
slowly varied). The Hamiltonian for this system in the 
rotating frame with U = cxp(— iujocr z t) and after making 
the optical rotating wave approximation is 

= n 12 {e l{ ^ ut+m2{a+af)) a + + h.c) 

+ f7 3 4(e^ t+J)34(a+Q+) V + + h.c) (12) 

where and a are the creation and annihilation opera- 
tors for the vibrational modes, cr_ = |g)(e|, er+ = |e)(g| 
and g z = \g)(g\ — |e)(e|, are Rabi frequencies of 
|e) — * \g) transition induced by the ith and jth lasers 
and the factor rjij = Skija Q is the Lamb-Dicke parameter 
(where a Q is the amplitude of the ground state of the trap 
potential and Skij is the wave vector difference between 
the ith and the jth beams). Hamiltonians comprising 
of any one of terms ^i l2 ^{<t>-^+ m ^a+a J ')) a+ + QJ . 

^34(e l(lyt+,,34(a+Q+)) cr + + h.c) have been used in the con- 
text of nonclassical state preparation of the vibrational 
mode |fil| and in a recent experiment such terms where 
switched on simultaneously p3| ]. These are generally 
implemented with a single pair of travelling wave laser 
beams. We require two pairs of laser beams to switch 
on both the Hamiltonian terms simultaneously. Simi- 
lar Hamiltonian terms, when used simultaneously and in 
conjunction with atomic decay, can be used to generate 
squeezed states of motion of the ion |l^j2^]. Here, how- 
ever, we do not want atomic spontaneous emission, and 
thus choose the |e) and |<7) to be hyperfine levels so that 
atomic decay can be neglected. In the Lamb-Dicke limit 
(riij « 1) the field can be expanded to the first order in 
rjij. Expanding thus, and transforming into the rotating 
frame with U = exp(— iva^at) we obtain 

ijM ^U^HU = \e)(g\(g a e l 'l'a + g b a t ) + h.c. 

= H a + H b (13) 

where we have made another rotating wave approxima- 
tion by omitting all the rapidly oscillating terms. In the 
above, g a — 7/12^12/2 and g b = 7734^34/2. Thus we can 
obtain the Hamiltonian H in an ion trap by the above 
methods in a rotating frame. From an experimental per- 
spective, instead of slowly varying the phase of one of the 
beams relative to the rest, one can alternatively keep the 
phase of all lasers same but set wli — ^L2 =w + ^ + (5 
where 8 « v. Then td which varies very slowly with 
time t will serve as the varying phase <fi. There is also an 
alternative way to implement the required Hamiltonian 
which decreases the number of laser beams to two. In 
this case, each laser drives a quadrupole transition such 
as the narrow S'1/2 — * -D5/2 transition in 40 Ca + [?] and 
the lasers would have to be oppositely detuned from the 
transition by v. 

We now describe our proposal to measure a Harmonic 
oscillator Berry phase in a vibrational mode of a trapped 
ion. It is assumed that the atom is cooled to the ground 
state of motion and in the internal state | g) . First we pre- 
pare a coherent state \a) of the vibrational mode. This 
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can be done by shifting the centre of the trap |L2| or 
other methods Jl!||l3| (our scheme is independent of the 
method used). Next we want to achieve the evolution 



lfl>l«>-lfl>l 



(14) 



by purely geometrically. To this end we switch on 
the Hamiltonian H by switching on appropriate lasers 
and viewing the dynamics from the rotating frame with 
U = exp[— i(ya^ a + u>Q<j z )t] (we will describe later what 
happens in the laboratory frame). The initial state can 
be written as 



k»l 9 > = — H 



i*,t> + !*;»■ (is) 



Next, the relative phase <j) of the lasers is varied 
very slowly and cyclically while the parameter r = 
tanh -1 gb/ g a is kept constant. If (f> is varied slowly 
enough so that adiabaticity holds, the state evolves as 

„-|a| 2 /2 n 

+ e^n(*)|*-( e ))), (16) 



where j n {4>) an d Pni'f) are geometric and dynamical 
phases respectively and e = re~ 1 ^ . After <j> has com- 
pleted an entire cycle (i.e. <j> = 2tt), we apply a state 
dependent phase shift \g) — > —\g) to the ionic state (by 
applying a 2n pulse Jl^l ) and this has to be done much 
faster than the evolution timescale of the system. This 
converts |\l/+(e)) — » |W~(e)) and vice versa. Now we vary 
the phase difference <j> again from 2ir to 4tt. The resulting 
state at the end of this second cycle is 



I*(4tt)> 



V2 

l*n(0» 



E 



a n e i(7n(47r)+/3+(27r)+/3-(27r)) 



(1^(0) 

(17) 



As /3+ = — f3~, the dynamical phase completely cancels. 
If we choose sinh 2 r = 1/4, we have 7„(47r) = — mr. Un- 
der such circumstances, the final state would be 



|*(4tt)) = -^=^E L W(K«» + l*n(0» 



V2 
= lff>|-«) 



(18) 



So the detection of the Berry phase now amounts to dis- 
tinguishing between the coherent states \a) and | — a). 
To this end, after switching off the adiabatic evolution, 
the entire state is given a negative displacement of —a. 
This reduces our problem to distinguishing between JO) 
(a vibrational state with no excitation) and | — 2a) [ p5[ . 
To accomplish this, the ionic internal states are allowed 
to interact with the vibrational mode by a Jaynes Cum- 
mings interaction |E5[ . In the case of no Berry phase, the 



probability of finding the ion in an excited state at any 
time t is P eo — 0, while, in the presence of a Berry phase 
the same probability is given by 



Jr e, — 2a — e 



^ (~2a)- 



n 



■ sm 



(19) 



where fl n +i is the Rabi frequency corresponding to an 
excitation number n+1. Note that the above method is 
not a perfect discrimination (it is impossible in principle 
as |0) and | — 2a) are nonorthogonal) . Incidentally, the 
squeezing in our proposal (r = 0.48) is smaller than the 
value r ~ 1.5 required for the four squeeze nonadiabatic 
geometric phase proposal fl7f| . If the original coherent 
state amplitude is \a\ ~ 1, we can ensure that this squeez- 
ing does not bring in contribution from the higher order 
terms in the expansion of Eq. ( |l2l) (From Ref . (2l) one can 
show that the probability of the state |3) in squeezing of 
|1) with r ~ 0.48 is already small enough to make the 
ratio of higher order terms to the first order term lower 
than 10~ 2 rj 2 ). One alternative to using a coherent state 
is to start with the superposition |0) + |1) (which should 
not be too difficult to prepare |2q] ) and convert it to the 
state |0) — |1) by following an identical procedure to that 
described above for coherent states. 

In the laboratory (nonrotating) frame, there will 
be an extra phase equal to ^ c de{^!^\S{e)^ {va^ a + 
uJoO' z )S(e)\^^) . For our choice of sinh 2 r = 1/4 this 
phase turns out to be 3vT where T is the time required 
to complete one cycle of the phase cf>. We can choose T 
in such a way that this phase becomes a multiple of 2ir. 
This will keep the magnitude of the phase unaltered from 
that in the rotating frame. 

Let us now examine the feasibility of our experiment 
with existing ion trap parameters. For adiabaticity, 
we require the time scale T of variation of the rela- 
tive phase <p to be greater than the dynamical timescale 
of the problem. The dynamical timescale is given by 
{Vtij^ij /2)~ 1 . We first consider experiments with 

ion QUI where rj - 0.2 and n/2ir - 500 kHz. 
The dynamical timescale, required to be smaller than T, 
is then about 0.33 x 10~ 5 s. T, on the other hand, is 
required to be lower than both the lifetime of the excited 
state (which can be up to 10 s p0[|) and the motional 
decoherence timescale (which is about 10 ~ 4 s |llj for 
9 Be + experiments) so that our assumption of neglect- 
ing all types of decoherence holds. If we set T ~ 10 -5 s 
(much less than both the decoherence time scales), then 
it is about thrice the dynamical timescale and the as- 
sumption of adiabaticity should hold (one could even try 
to ensure an order of magnitude difference between the 
dynamical timescale and by increasing the laser power a 
bit more than three folds). On the other hand, motional 
decoherence, which is an order of magnitude slower than 
T, can be neglected. In 40 Ca + experiments Q, the 
dynamical timescale is 10 -5 s whereas the motional de- 
coherence time scale is 10 -3 s (internal state decoherence 
being 1 s). For these traps, we can choose T ~ 10 -4 s 



3 a lb 
9 Be 
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and perfectly satisfy both adiabaticity and the neglect of 
decoherence. We also have to set 3z/T — m2ir for the 
extra phase in the nonrotating frame to vanish, where m 
is an integer. For standard traps, v ~ 10 MHz, jl9| and 
hence vT ~ 10 4 . We can easily choose a number of such 
a large order to be a multiple of 2tt. 

It is interesting to point out that if our scheme is per- 
formed starting with the ion in the excited state, the 
evolution would be \e)\a) — > — |e)| — a). Therefore, the 
initial internal state of the ion makes no difference to 
our scheme. Using a similar scheme we can also create a 
Schrodinger cat state of the ion. One initially has to pre- 
pare the ion in a superposition of \g) and some other state 
\r) which is completely decoupled from the evolution due 
to our Hamiltonian. Then the evolution proceeds as 



(\9) + \r))\a) 



x) + \r)\a) 



(20) 



The state in the above equation is a Schrodinger cat state 
of the ion @. 

We have shown how to observe the Berry phase of 
a simple harmonic oscillator using a trapped ion. We 
have described how to reverse the amplitude of a coher- 
ent state by purely geometric means (i.e. by using only 
the Berry phase) . An important advantage of the exper- 
iment is that it needs only a single trapped ion. This is a 
simple requirement in comparison to the technology that 
already exists such as the ability to entangle four ions in 
a trap ^3|. Moreover, we have shown that the existing 
ion trap parameters are well in range of those required for 
implementing our proposal and for our analysis to remain 
valid. We have also indicated how to use our scheme to 
create a Schrodinger cat state of the ion. Further inter- 
esting extensions to the geometric approach to multiple 
modes are also possible and will be investigated in the 
future. 
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